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Abstract 

In previous publications dressed coordinates and dressed states has been intro- 
duced in linear systems. Specifically, a system composed by a harmonic oscillator 
interacting linearly with an infinity set of other oscillators has been treated. In this 
paper we show how to generalize such dressed coordinates and dressed states to a 
nonlinear version of this system. Also we clarify some misunderstandings about the 
concept of dressed coordinates. Indeed now we prefer to call them renormalized 
coordinates to emphasize the analogy with the renormalized fields in quantum field 
theory. 
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1 Introduction 



In recent publications it has been introduced the concept of dressed coordinates and 
dressed states[l, 2, 3, 4], in the context of a harmonic oscillator (the atom) interacting 
linearly with a massless scalar field. This allowed the construction of dressed states, with 
the physically correct property of stability of the oscillator ground state in the absence of 
field quanta (the quantum vacuum). For a review see the next section. Also this formalism 
allowed a nonperturbative treatment of the oscillator radiation process. When applied to 
a confined atom (approximated by the oscillator) in a spherical cavity of sufficiently small 
diameter the method accounts for the experimentally observed inhibition of the decaying 
processes [5, 6]. 

In Ref. [7] an attempt to construct dressed coordinates and dressed states for a 
nonlinear system has been done. However the approach used there was more intuitive 
than formal. The purpose of this paper is to develop a formal method to construct 
dressed coordinates in nonlinear systems. We will do this by a perturbative expansion in 
the nonlinear coupling constant. To be specific we consider the model with Hamiltonian 
given by, 



1 1 

H = 2 ( p o + u b<1o) + « E {Pk + "fal ~ 2c k q k q ) 

k=l 

N N 
r=0 r=0 

where sums over repeated indices and the limit iV — > oo are understood. In Eq. (1) the 
bare frequency of the oscillator, cu B , is related to the physical frequency, u , by [8, 9], 



N „2 



2 > t—v C 



o + ES- ( 2 ) 



k=i u k 



The coefficients T^) and TZ^ paTe are chosen in such a way that the Hamiltonian given 
by Eq. (1) is positive defined. 

In Ref. [7] the quartic nonlinear model, a r = 0, has been treated. Here also we will 
be mainly interested in the quartic nonlinear model to compare with the early treatment. 



2 



The sextic nonlinear interaction will be considered mainly because for some values of A r 
and a r it is possible to find the exact solution for the ground state and, as explained at the 
end of section III, this fact will permit an exact construction of the dressed coordinates. 
Then, this sextic model will allow us to test the validity of the strategy developed to 
obtain the dressed coordinates in a general nonlinear system. 
Along this paper we use natural units c — h — 1. 

2 Defining dressed (renormalized) coordinates and 
dressed states 

The purpose od this section is twofold. First, to make this paper self contained we review 
what was called dressed coordinates and dressed states in the Refs. [1, 2, 3]. Second, we 
clarify some misunderstandings about the concept of dressed coordinates, as we explain 
below. Indeed now we prefer to call them renormalized coordinates to emphasize that 
these coordinates are analogous to the renormalized fields in quantum field theory. 

To understand the necessity of introducing dressed coordinates and dressed states let 
us consider the following problem. Imagine that the oscillator with index zero in Eq. (1) 
represents an atom and the other oscillators represent the modes of the electromagnetic 
field. If there are no interaction among them the free Hamiltonian (cj, = X r = a r = 0) 
has the following eigenfunctions, 



N 

= nvvM, (3) 

where \q) = |?o, 9i, 9jv) an d Vv^m) ^ s ^ ne eigenfunction of a harmonic oscillator of 
frequency u^, 



nom...nN 



(q) 



(g|n ,ni, ...,n N ) 
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The physical meaning of il> no n 1 ...n N (<l) in this case is clear, it represents the atom in its uq- 
th excited level and n k photons of frequencies lo^. Now, consider the state of no photons, 
VW)-o(<z) : the excited atom in the quantum vacuum. We know from experience that 
any excited level of the atom is unstable, that is, the state VVi o...o(9) is unstable. The 
explanation of this fact is that the atom is not isolated from interacting with the quantum 
electromagnetic field, or in other words it is rendered unstable by interacting with the 
quantum vacuum. This interaction in our model is given by the linear and nonlinear 
couplings of go with q k . Obviously, when we take into account these interactions any state 
of the type VW)-o(9) wm be unstable, since these states are not eigenfunctions of the total 
interacting Hamiltonian. But, there is a problem, the state ^oo.-oC?)) that represents the 
atom in its ground state and no photons, is also unstable contradicting the experimental 
fact of the stability of the ground state in the absence of photons. What is wrong in all 
this? The first thing that cames in our mind is to think that the model given by Eq. (1) 
is wrong. Certainly, we know that the correct theory to describe this physical system is 
quantum electrodynamics. On the other hand such a description of this system could be 
extremely complicated. If we aim to maintain the model as simple as possible and still 
insist in describing it by the Hamiltonian given in Eq. (1) what we can do in order to 
take into account the stability of the ground state? The answer lies in the spirit of the 
renormalization program in quantum fiel theory: the coordinates g M that appear in the 
Hamiltonian given by Eq. (1) are not the physical ones, they are bare coordinates. We 
introduce renormalized coordinates, q' and q' k , respectively for the dressed atom and the 
dressed photons. We define these coordinates as the physically meaningful ones. These 
renormalized coordinates were called in prededing works as dressed coordinates [1, 2, 3], 
for this reason, from now on we will take these denominations as synonymous. In terms 
of these dressed coordinates we define dressed states as 



/i=0 

where the subscript d means dressed state, \q f ) = \q' ,q[, ■ ■■,q'iy) and ipn^Qu) is given by 
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(6) 



The dressed states given by Eq. (5) are defined as the physically measurable states and 
describe in general, the physical atom in the no-th excited level and rik physical photons of 
frequency uj^. Obviously, in the limit in which the coupling constants Ck, A r and a r vanish 
the renormalized coordinates q' approach the bare coordinates q^. Now, in order to relate 
the bare and dressed coordinates we have to use the physical requirement of stability of 
the dressed ground state. The dressed ground state will be stable only and only if it is 
defined as the eigenfunction of the interacting Hamiltonian given in Eq. (1). Also the 
dressed ground state must be the one of minimum energy, that is, it must be defined 
as being identical (or proportional) to the ground state eigenfunction of the interacting 
Hamiltonian. From this definition, one can construct the dressed coordinates in terms of 
the bare ones. 

Firstly we explicitly construct the dressed coordinates for the linear model obtained 
from Eq. (1) by setting \ r = a r = 0, 



Although the task of constructing dressed coordinates in linear systems has been done in 
preceding works, we repeat here the calculation in order to make this paper self contained. 
In the next section we will consider the nonlinear case. As is well know, the Hamilto- 
nian (7) can be diagonalized by means of the introduction of normal coordinates Q r and 
momenta P r , defined as 





N N 



9, = E«;a, ?v = E^' ^ = (o,*o, k = i,2,...,N, 



(8) 



r=0 r=0 



where {t r } is an orthogonal matrix whose elements are given by [10, 11] 
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(9) 



In normal coordinates the Hamiltonian (7) reads 
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1 N 



H linear = -J2(P? + n 2 r Q 2 r ), (10) 

A r=0 

where the fi r 's are the normal frequencies, corresponding to the collective modes and 
given as solutions of [10, 11], 

N 2r>2 

-I - ^ 2 = E ■ (») 

fc=l ^fel^fc _ "rJ 

The eigenf unctions of the Hamiltonian given in Eq. (10) are given by 



(Q) = {Q\n ,n 1 ,...,n N ) l 

N 



= II^(Qr), (12) 



r=0 



where the subscript c means collective state, \Q) = \Qo,Qi, ■■■,Qn) and (p nr (Qr) are the 
wave functions corresponding to one dimensional harmonic oscillators of frequencies Q r , 

^^^"EMl^. (13) 

Now, the dressed coordinates are defineed requiring t/>oo...o(<?') ^ (j>oo...o(Q), since in this 
way we guarantee that tpoo...o(q') is the ground state of Hu near . Then, from Eqs. (12) and 
(5), we have 



e 4E M N = 0^«) 2 Qc e -^Er=0 n rQ 2 r ? ( 14 ) 



from which we obtain 



^ = E J^Wr , (15) 

as can be seen by direct substitution in Eq. (14) and using the orthonormality properties 
of the {t^} matrix. The above definition guarantees the stability of the dressed ground 
state, however, since the other dressed states are not energy eigenfunctions, they will 
not remain stable. For example the first excited dressed state, whose eigenfunction is 
V'io...o(<Z / )> w iH decay to the ground state ?/>oo...o (<?')• 
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We have to remark here that the dressed coordinates here introduced are not simply 
a change of variables, they are new coordinates in its own right and are introduced by 
physical consistence requirement of the model. Also we have to stress that our dressed 
coordinates are not the same as the ones employed in other references, as for example 
in [12] and references therein, where the authors called dressed coordinates the collective 
normal ones. Also our dressed states are very different from the ones defined in Refs. 
[13, 14, 15, 16], where the authors called dressed states the states obtained by diagonalizing 
a truncated finite matrix representation of the Hamiltonian. 

Before leaving this section it will be useful to establish the relation between 
i>nom...n N (q') = ■-,n N ) d and {Q\n Q ,n u ...,n N ) d . For this end we write 



i{no,n 1 ,...,n N \m ,m 1 ,...,m N )d = J dq' d (n , n 1: n N \q'){q'\m , m u m N ) 

dq' 



J dQ 



dQ 



i(n ,ni, ...,n N \q')(q'\m ,mi, ...,m N ) d 



= J dQ d (n ,ni, ...,n N \Q)(Q\m ,m 1 , ...,m N ) d , (16) 

where dq' = Y\^ =0 dq' dQ = Y\^ =0 dQ r and \dq'/dQ\ is the Jacobian associated to the 
transformation q'^ — > Q r . From Eq. (16) we get 



{Q\no,ni, ...,n N ) d 



dq' 



dQ 



1/2 



(q' {710,71!, ...,n N ) d ■ 



(17) 



Taking n = n\ 



n N = in Eq. (17) and using ■0oo...o(?') oc </>oo...o(Q) we get 



|0,0,...,0) d cx J dQ 



dq' 



dQ 



1/2 



|Q)(Q|0,0,...,0) c . 



(18) 



In the linear case, we easily get, from Eq. (15), \dq'/dQ\ = n^=o \J^r/^n and using this 
result in Eq. (18) we obtain 



|0,0,..,0> d oc |0,0,...,0) c . (19) 

For a nonlinear system, certainly a relation of the type given by Eq. (19) will not hold. 

In next section we construct dressed coordinates in the nonlinear model described by 
the Hamiltonian given in Eq. (1). 
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3 Constructing renormalized coordinates in a non- 
linear model 



Now we are ready to construct dressed coordinates and dressed states in the nonlinear 
model given by Eq. (1). For this purpose we have to find, firstly, the eigenfunctions of 
this Hamiltonian, in particular its ground state eigenf unction. 

In order to maintain things as simple as possible and to compare with the the preceding 
treatment given in Ref. [7], we consider the nonlinear quartic interaction obtained from 
the model given described in Eq. (1) by setting a r = 0. Following Ref. [7] we make the 
simplest choice for the coefficients as 

T\xvpo = tptjstpte . (20) 

Substituting Eqs. (8) and (20) in Eq. (1) we get 

1 N 

H=-J2(P? + n 2 r Q 2 r + 2\ r Qt) , (21) 

that is, we obtain a system of uncoupled quartic anharmonic oscillators. In Eq. (21) we 
can notice that A r has dimension of [frequency] 3 . Then we write A r = \Q 3 , where A is 
a dimensionless coupling constant. The eigenfunctions of the Hamiltonian given by Eq. 
(21) can be written as 



<j>nom...n N {Q; A) = (Q\n ,n 1 ,...,n N ;X) c 

N 

= Y[<Pn r (Q r ;X) , 

r=0 



(22) 



where 4> nr (Qr', A) are eigenfunctions of (P r 2 + Q 2 .Q 2 + 2\Q 3 Qj) /2 and can be written for- 
mally as (see Appendix) 



*>n r (Qr; A) 



7T 



Hn r ( y/QrQr ) 



i=i 



g 2 



(23) 



where G$(yJTL r Q r ) are linear combinations of Hermite polynomials. The corresponding 
eigenvalues for the Hamiltonian given in Eq. (21) are given by, 
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N 



E, 



non±...riN 



(24) 



r=0 



where E nr (X) are the eigenvalues corresponding to the eigenstates given in Eq. (23), 



1 oo 

E nr (X) = (n r + -)n r + J2\ l E£, 
z i=i 



(25) 



with the obtained by using standard perturbation theory (see Appendix). 

Taking no = n\ = ... = = in Eq. (22) we get for the ground state eigenfunction 
of the total system, 



N /ft 

0oo...o(g;A) = n — 

r=0 V 7F 



2 



(26) 



To properly define (see comments below) the dressed coordinates it is convenient to write 
the above equation as, 



W /n \ 1/4 T oo oo / 

o(0;A) = n(-) i+E^ + E^F^)-^ ) 

r=0 ^ 71 ' L 1=1 1=1 ^ 



^00.. 



g 2 v ' 



N 

oc n 

r=0 



n=0 /o/l.../n=l 



e-^, (27) 



x (G$°\ v^Q r ) - } (0) J G^O)...^^) 
where in the second line we factored the term 1 + £~i \ l G$\Q) and used (1 + a;)" 1 = 

££*(-i) B * B . 

The physically measurable states, the dressed states, are defined by Eqs. (5) and 
(6). Hence, the dressed coordinates q'^ will be defined in such a way that the dressed 
ground state equals (or is proportional) to the ground state of the nonlinear interacting 
Hamiltonian given in Eq. (27). That is, we define the dressed coordinates imposing the 
condition ipoo...o(q') oc 0oo...o( < 3; A) which by using Eqs. (5), (6) and (27) can be written as 



1 V^-W / / \2 



r=0 L ra=0 l h...l n =l 

x (G$°\JsKQ r ) - gS o) (0)) G^(0)...G^\0) 



e 2 



SJ fQ 2 r . (28) 
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Now, we write a perturbative expansion in A for q' , 



N rrr 



r=0 V U ^ 



i oo 

Q r + —£\*FW(y/n r Q r 
V "r i = i 



(29) 



Replacing Eq. (29) in Eq. (28) and using the the orthonormality of the matrix {t^} we 
get 



cxp 



OO 1 oo 

£\ l y/n r Q r FP(y/n r Q r )- - Yl \ l+m F^ l \^ r Q r )F^ m \^ r Q r ) 

1=1 ' ^ l,m=l 



= 1 + El" 1 )" E A ,0+,1+ - +,n (G? o) (v^rQr) - G? o) (0)) G^\0)...G { ^\0) . (30) 

n=0 i li...! n =l ^ ' ' 

Expanding the exponential in the left hand side of Eq. (30) and identifying equal powers 
of A we can obtain all the F^{\fW r Q). The general expression is very complicated, here 
we only write the first three terms: 

^r (1) (er) = -^(Gf ) (er)-Gf ) (0)) , (31) 
Sr 

F?\tr) = ~\ \d$\ir) ~ Gf(0) + G$-\0)t r FV(t r ) + hi - £) (F«(£ r )) 2 ] (32) 



and 



G?\tr) ~ Gf (0) - (<#>(£.) - Gj 2) (0)) G?>(0) 

3! 



+ trFUfr) (g$\o) - (G«(0)) 2 ) + I (C 2 - 3) 5 
+(l-&)FV(Z r )FW(Z r )] , 



(33) 



where £ r = y/Q.Q r . From the Appendix, using Eqs. (68-71) in Eqs. (31) and (32) we get 
respectively, 



(34) 



and 
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^r) = -^(y& + 14£ + ^ • (35) 



Replacing above equations in Eq. (29) we obtain at order A 2 , 



N 



r=0 



(36) 



where we have introduced the dimensionless dressed coordinate = y/UJ^q'. 

Before leaving this section we would like to comment why we factored the term 
1 + ^ 1 Gq\o) in Eq. (27). Notice that we define the dressed coordinates q'^ by means 
of the proportionality 1P00...0W) 4>oo...o{Q] A). To promote this proportionality into an 
equality we have to take care in defining a well behaved transformation between dressed 
and collective coordinates, for example it would be undesirable any singular transforma- 
tion. To see how the above undesirable situation occurs, define the dressed coordinates 
through Eq. (26) without the factorization of the term 1 + J2tli ^ 1 Gq\o). It is easy to 
show that in such a case the transformation between q'^ and Q r is singular. For example 
we obtain for F^(^ r ) and F r 2 (£ r ) 

F?\Zr) = -±-<$ ) {Zr), (37) 



(38) 



2 

Since G [ L> (^ r ) and G^'fa) are n °t homogeneus functions of £ r [see Appendix, Eqs. (68) 
and (69)] then Eqs. (37) and (38) are singular in £ r = 0. Consequently, the dressed 
coordinates defined through this prescription are not well defined. To understand what 
is happening and how to remedy this problem, note that this singularity means that 
€ r Fr\€r) is not homegeneus in £ r . But the effect of this nonhomogeneous term on the 
wave function [that contains terms of the type e _A ^ rFr( } see Eq. (30)] is just equal 
to a factorization term. Then to remedy the situation we have to make a convenient 
factorization in ■0oo..o(?') or, equivalently, in 0oo..o(Q) before promoting the proportionality 
into an equality. That our choice, the factorization of 1 + ^ l Go\o) in Eq. (27), is 
the correct one is supported by the fact that we obtain well behaved dressed coordinates. 
To further support our choice, we consider a system in which we can solve exactly for 
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the ground state, allowing us to obtain exact dressed coordinates. Comparing these exact 
dressed coordinates and the perturbative ones we get the same answer. The model in 
mention is the one whose Hamiltonian is given by Eq. (1) with coupling constants defined 
as 

AO 3 A 2 fi 4 

\ r T{ r ) r j j .r j ip( r ) r .r j.r j_r j.r yr +r ("lO} 

A r I p J upa ~ ^ _ Jj^3/2 V " V r ^fivpare ~ _ 3_\)2 P " P u T e ' ^ ' 

Using the above expression in Eq. (1) we get a system of uncoupled sextic anharmonic 
oscillators, 

1 N ( 9X0 3 X 2 4 \ 

H S (* + ^ + <^w Q * + t^w q ") ■ (40 » 

By direct substitution it is easy to show that the above Hamiltonian have as ground state 
eigenfunction the following solution [17] 

0(Q; A) = Me-^ N r=^ +x ^) , (41) 
where M is a normalization constant, 

ft = vra < 42 > 

and the corresponding ground state energy is given by 

£(A)=f>. (43) 

r=0 

Now, the dressed coordinates can be defined by 



from which we obtain 



' 1 k v 2 



(45) 



^ r ? ^ r (v / T^3A + 2(l-3A) 
Note that at order A both the quartic and sextic anharmonic Hamiltonians, given respec- 
tively by Eqs. (21) and (40), are equivalent. Then, if our strategy to define perturbatively 
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the dressed coordinates is correct, at order A Eq. (36) must be identical to Eq. (45). Ex- 
panding Eq. (45) at order A we can see that it is indeed the case. Then we conclude that 
our strategy for defining perturbatively the dressed coordinates is the correct one. 

4 The decay process of the first excited state 

In Ref. [7] the probability of the particle oscillator to remain in the first excited state has 
been computed at first order for the nonlinear quartic interaction. However, as we have 
already mentioned, the approach used there was more intuitive than formal. In order 
to see in what extent such calculation is correct, in this section we compute the same 
quantity by using the formalism presented in the last section. To maintain the reasoning 
as general as possible we present the steps necessary to compute the probability amplitude 
associated with the most general transition, 

AZT^nTit) = d(m ,m 1 ,...,m N \e-' lHt \n ,n 1 ,...,n N ) d , (46) 

that is, we prepare our system initially at time t = in the dressed state \no,iii, ■ ■■,n^) d , 
then we ask what is the probability amplitude of finding, in a measurement performed at 
time t, the dressed state \m ,m 1 , ...,m N ) d . Introducing a complete set of eigenstates of 
the total Hamiltonian H, given by Eq. (22), in Eq. (46) we find 

oo 

d{m ,mi, ...,m N \e~ lHt \lo,h, —,In',^)cc{IoA, —,lN',M n o,n-i, — ,™Jv)d 

IoIi—In=0 
oo 

E ttil^fti^W^'^. ( 47 ) 

IoIi—In=0 

fti(A) = JdQ c (l ,l 1 ,...,l N ;MQ)(Q\no,n 1 ,...,n N ) d 

1/2 

<f)i oh ...i N {Q; (</) , (48) 

where in the second line we have used Eq. (17). 



Am m 1 ...m N (j\ 
• Al n a ni...n N \ L ) 



with 
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From Eq. (29) we get easily the Jacobian \dq' /8Q\, 



dq' 



dQ 



N 

n 

r,fi=0 



CO 



(49) 



Now we evaluate, at first order in A, the probability amplitude to the particle oscillator 
remain at time t in the first excited state if it has been prepared in that state at time 
t — 0. This quantity is obtained taking n = m = 1 and n k = m k = in Eq. (47). 
Notice that to compute -4-i§;;;§(t) at first order in A we need to find /jv (A), defined in 
Eq. (48), at order A. Replacing Eq. (34) in Eq. (49) we get, 



dq' 



dQ 



N n 

v/u,r=0 U t* 



1/4 



o\ N , s 

1 + 32 E [QHo(^sQs) + H 2 (^Q s Q s j 



+ 0{\ 2 ) . (50) 



At order A, from Eq. (22), we have for (/>i i 1 ...i N (Q, A), 



JV 



N 



<t>l h...l N {Q,X) = HMQr) + A]T 
r=0 r=0 

where the 4>i r {Qr) are given by Eq. (13), 



(51) 



G^UttrQr) = ai r H lr _ 4 (Jtt r Q r ) + bi r Hi r - 2 (\/tt r Q r ) + c lr H lr+2 (J tt r Q 



+di r Hi r+4: (J Q r Q r ) , 



(52) 



and ai r , bi r , c\ r and d\ r are given in Appendix, Eq. (70). Using Eq. (28) and Eq. (29) 
we have for ^io...o(</)> 



N 

n 

N 



1/4 



= n 



7T 



7T 



glVV^Ogo) e -^E; = o^(^) 2 

\/2 



1/4 



H 1 (y/n r Q r ) + 2XF^\y/n r Q r ) 

r=o v 2 



at , N 

- xh^qmy.^sQs FW{yJn s Q s ) e-^.=»^ + o(a 2 ) . (53) 



s=0 
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Replacing Eqs. (50), (51) and (53) in Eq. (48) we obtain after a long, but straightforward, 
calculation 



N - v/fi N 

Tlto lN W = £^<UlKo + 9^-A E*5*wIRo 

r=0 s^r r=0 s^r 

+3^AE^U^ 2 II Ko + 0(\ 2 ). (54) 

r^s Uy^r,s 

Replacing Eq. (54) in Eq. (47) and using from the Appendix Ei r (X) ^§f2 r +^Afi r we 
obtain for *4.}o;;;o(t), which we denote as foo(t; A), 

fait; A) = E^,nr J2(t r ) 2 (l - ^-i\tQ r ) + 0(\ 2 ) . (55) 

From the above equation we get the probability to the particle oscillator remain in the 
first excited level, 

1/ooC*; A)| 2 = |/oo(0| 2 + ^l|l/ooW| 2 + 0{\ 2 ) , (56) 

where 

N 

/«>(*) = E(^) 2 e-^. (57) 

r=0 

Equation (56) is the same as the one obtained in Ref. [7]. We obtained the same 
result because at order A the square of T^'q 1n is given only by the square of the first 
term in Eq. (54), that does not depend on A. Then, the effects of the nonlinearities, at 
this order, cames only from the corrections to the energies, as was assumed in Ref. [7]. 



5 Conclusions 

In this paper, after clarifying what we understand by dressed coordinates and dressed 
states, we have developed a formal method to construct perturbatively dressed coordinates 
in nonlinear systems. Although we restricted our calculations to a very special quartic 
interacting term, we have pointed out the necessity of factoring a term in order to avoid 
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artificial singularities which otherwise would appear if we do not make such factorization. 
That this factorization is the correct one has been checked by using an exactly solvable 
sextic interacting model. Then, in more general nonlinear systems, one can follow the 
same procedure to construct the dressed coordinates. 

At the end of section II we remarked that for nonlinear systems, in the number rep- 
resentation, the dressed ground state is not equivalent to the ground state of the total 
system, see Eq. (18). This fact must not be seen as in contradiction with our defini- 
tion of dressed coordinates, since we have defined them by requiring the equivalence of 
the dressed ground state in dressed coordinates representation and the ground state of 
the system in normal coordinates representation. We can understand the mentioned non 
equivalence, by noting that although the dressed ground state is an eigenstate of the 
dressed number operators (associated with the dressed coordinates) the ground state of 
the system, in general, is not an eigenstate of the collective (normal) number operators. 
For example, in the quartic nonlinear case one can easily verify that the ground state 
(and also, the other eigenstates) is not an eigenstate of the collective number operators, 
but a linear superposition of eigenstates of these operators [see Appendix, Eq. (62)]. 

Finally, we considered the computation of the probability of the particle oscillator to 
remain excited in its first excited level, and showed that the result coincides with the 
result obtained in Ref . [7] . Then, one of the conclusions of Ref . [7] remains: the effect 
of the nonlinear quartic terms is the enhancement of the decay of the particle oscillator 
from its first excited level to the ground state. This fact, can be easily seen from Eq. 
(56) by noticing that |/oo(£)| 2 (the probability to the particle oscillator remain in the first 
excited state in the absence of nonlinear interactions) in free space, is a decreasing (almost 
exponentially) function of time. 
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Appendix: The perturbed eigenfunctions and eigen- 
values 



It is easy to see that the eigenfunctions of the quartic anharmonic oscillator can be written 
formally as those given in Eq. (23). We have to notice only that any wave function can 
be expanded in the basis <p n (Q) (we omit here the index r), given by the eigenvalues of 
the linear part of the Hamiltonian. And since <fi n {Q) are given by exp(— QQ 2 /2) times 
a Hermite polynomial of degree n, we see that an expression of the type given in Eq. 
(23) follows. In what follows we compute G^(VttQ) and G^(VUQ) by using standard 
perturbation theory. 

At second order in standard perturbation theory the eigenfunctions and eigenvalues 
of a Hamiltonian H = H + XV are given respectively by, 



k ^ n ^n-^k \ K¥n [&n ~ E k ){E n ~ &l) k+n (A, ~ H) ) 

(58) 

and 



w, 



12 



E n {\) =E n + XV nn + A 2 £ -f^r + 0(\ 3 ) , (59) 

k+n Ai ^ k 

where 

V kn = (k\V\n) (60) 

and \n) and E n are respectively eigenfunctions and eigenvalues of the unperturbed Hamil- 
tonian H . 

For the anharmonic oscillator with V = il 3 Q 4 we obtain easily 



kn 



j[y/hS k , n -A + 2(2n- l)Jk 2 Sk,n-2 + 3(2n 2 + 2n + l)S k:T 



(61) 
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where k n = (k + l)(k + 2)...(k + n). Replacing Eq. (61) in Eqs. (58) and (59) we obtain 
respectively, 

|n,A) = \n) + X(a' n \n- 4) +b' n \n-2) + c' n \n + 2) + d! n \n + 4)) 
+ A 2 (e' n \n - 8) + f n \n - 6) + g' n \n - 4) + ti n \n - 2) 

+<Jn + 2) + u' n \n + 4) + <|n + 6) + iwjjra + 8)) + £>(A 3 ) (62) 



and 



where 



E n (X) = (n + + A£« + A 2 £f + 0(A 3 ) , (63) 



< = -^f^, < = -^; (64) 



512 W- ,n lg2 

= 1(2^-9^ + 7)^^-4)4, ^ = -^(2n 3 + 93n 2 - 107n + 66)^(71 - 2) 2 , 

< = -l(2n 3 -123n 2 - 359n-300)v^, < = l(2n 2 + 13n + 18)^ , 

(6n + 17) 1 
U " = 192 — Wn= 512 v ^' ^ ^ 

= |j(2n 2 + 2n + l)ft (66) 

and 

= --(34n 3 + 51n 2 + 59n + 21)fi . (67) 
8 

Writing Eq. (62), in coordinate representation, in the form given in Eq. (23) we get for 
G^(VUQ) and G^(V^Q) respectively, 
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= On#n-4(0 + 6n^n-2(0 + Cn#n +2 (0 + d n H n+4 (£) (68) 

and 

Gf(0 = e n H n _ 8 (C) + f n H n _ 6 (0+9„H n _,(0 + h n H n _ 2 (0 

+t n H n+2 (£) + U n # n+4 (0 + U n -ffn+6(0 + W n H n+8 (£) , (69) 

where £ = VTlQ, 



v /2"- 2 (n-2)! ' 

r„ = , '" d B = d>n (70) 

' ' ^2"+ 4 (n + 4)! 



and 



a 'n 


^2 n ~\n 


-4)! 


d 




^2»+ 2 (ra 


+ 2)! 








-8)! 


g' n 




^2 n ~ 4 (n 


-4)! 






v /2™+ 2 (n 


+ 2)! 






v /2"+ 6 (n 


+ 6)! 



f 

J n 

n 



9n / ? 



t = V 

Vfi t — — - ■ — , u r 



^2 n - 6 (n 


-6)! 


K 




^2 n ~ 2 (n 


-2)! 


< 




v /2«+ 4 (n 


+ 4)! 


< 




v /2"+ 8 (n 


+ 8)! 



W n = , . (71) 

v /2"+ D (n + 6j! 
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